EPJ manuscript No. 

(will be inserted by the editor) 



Dark lattice solitons in one-dimensional waveguide arrays 
with defocusing saturable nonlinearity and alternating 
couplings 

Andrey Kanshu 1 , Christian E. Riiter 1 , Detlef Kip 1 , Jesus Cuevas 2 , and Panayotis G. Kevrekidis 3 

1 Faculty of Electrical Engineering, Helmut Schmidt University, 22043 Hamburg, Germany 

2 Grupo de Ffsica No Lineal, Departamento de Ffsica Aplicada I, Universidad de Sevilla. Escuela Politecnica Superior, C/ 

Virgen de Africa, 7, 41011 Sevilla, Spain 

3 Department of Mathematics and Statistics, University of Massachusetts, Amherst MA 01003-4515, USA 

May 4, 2012 

Abstract. In the present work, we examine "binary" waveguide arrays, where the coupling between adjacent 
sites alternates between two distinct values C\ and C2 and a saturable nonlinearity is present on each site. 
Motivated by experimental investigations of this type of system in fabricated LiNb03 arrays, we proceed 
to analyze the nonlinear wave excitations arising in the self-defocusing nonlinear regime, examining, in 
particular, dark solitons and bubbles. We find that such solutions may, in fact, possess a reasonably wide, 
experimentally relevant parametric interval of stability, while they may also feature both prototypical 
types of instabilities, namely exponential and oscillatory ones, for the same configuration. The dynamical 
manifestation of the instabilities is also examined through direct numerical simulations. 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

Dark solitary waves are ubiquitous nonlinear excitations 
of dispersive wave models with a so-called self-defocusing 
nonlinearity. Perhaps the prototypical model where they 
arise is the nonlinear Schrodinger equation with a defocus- 
ing nonlinearity pQ, as it is referred to in nonlinear optics 
(since it induces the spreading of optical beams). These 
excitations consist of a density dip (i.e., a dark notch) ac- 
companied by a phase jump across the density minimum. 
In addition to their observation in nonlinear optics exten- 
sively summarized in [UH], numerous experiments have 
demonstrated the emergence of these nonlinear states in 
the atomic physics of Bose-Einstein condensates, as has re- 
cently been comprehensively reviewed in [3] . Interestingly, 
the physical relevance of dark solitons are not limited 
to these areas but rather also extend to parametrically- 
driven shallow liquids [4] , discrete mechanical systems [5] , 
electrical lattices with nonlinear capacitors [6], thin mag- 
netic films [7, and dissipative variants thereof in complex 
plasmas [8], among others. 

In recent years, there has been a considerable interest 
in the realization of such excitations in systems bearing 
lattice potentials, to appreciate the effects of discreteness 
on the existence and stability, as well as on the dynamics 
and mobility of such excitations; see e.g. the prototypical 
theoretical study of [9] and even the earlier work of [TO] . 
Experimental developments have enabled the realization 



of such excitations, especially so in arrays of optical waveg- 
uides either in the context of AlGaAs in the anomalous 
diffraction regime with the cubic nonlinearity due to the 
Kerr effect [IT], or in defocusing lithium niobate waveg- 
uide arrays, which exhibit a different type of nonlinearity, 
namely a saturable, defocusing one due to the photovoltaic 
effect [12] . These experimental investigations, in turn, led 
to comparative studies between the features of the dark 
solitons in these two discrete models [13,14 , as well as to 
the examination of states consisting of multiple dark soli- 
tary waves [T5| [T6] . These findings have been summarized 
in a recent book [17] . Furthermore, in these waveguide ar- 
ray systems, not only have dark solitons been identified in 
higher gaps (of the associated periodic potentials) [18] , but 
they have also been found to arise as members of multi- 
component soliton complexes such as dark-bright solitary 
waves [T9] . 

In the present work, we consider a nontrivial variant 
of the above nonlinear dynamical lattices. In particular, 
we focus on the setting of "heterogeneous" lattices and 
more specifically in the context of lattices with alternat- 
ing couplings between two distinct values C\ and C2. The 
existence and stability of bright solitons in binary lattices 
has been investigated both theoretically [20,21 and ex- 
perimentally [22,23]. It should be added, however, that 
binary lattices can be constructed in different ways, such 
as e.g. by waveguides of the same separation between the 
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channels, but alternating in widths as in [22], or by ones 
of the same width but alternating in separation as in [23] . 
Here, in line with the latter approach, we fabricate one- 
dimensional arrays of such lattices in a nominally non- 
doped LiNbC>3 substrate by Ti in-diffusion. In this case 
the underlying linear spectrum possesses two pass bands, 
with a mini-gap between them (as opposed to the single 
band of the homogeneous waveguide array). In the gaps 
of the linear spectrum, we will seek nonlinear wave solu- 
tions, such as dark solitons and bubbles in what follows. 
We find that both types of obtained nonlinear waves may 
possess both regions of stability as well as regions of both 
types of instability, exponential as well as oscillatory, as 
the propagation constant parameter is varied. This is con- 
trary to what is the case for the homogeneous (coupling) 
counterpart of the model, where only oscillatory (or only 
exponential) instabilities may be observed for a particular 
type of configuration. When the waveforms our found to 
be unstable, direct numerical simulations of the system 
are used to explore the evolution of the instability, giving 
rise to the breakup of the stationary structure into moving 
states or to its breathing dynamics. 

Our presentation is structured as follows. In section II, 
we present the experimental setup and some prototypical 
experimental results. These, in turn, motivate the theoret- 
ical model presentation and discussion of section III and 
the systematic numerical investigation of section IV. Fi- 
nally, in section V, we summarize our findings and propose 
some potential directions for future study. 



2 Experimental Motivation 

For our experiments, we fabricate a one-dimensional waveg- 
uide array in a nominally non-doped LiNbOa substrate. 
Using standard photolithographic techniques, we form, on 
the substrate's surface, an array with 250 Ti stripes of 
equal width W = 4.0 jam and alternating distances d\ = 
2.5 /im and =4.5 /am between them. By in-diffusion of 
Ti stripes at T = 1000 °C and following polishing, we ob- 
tain samples with dimensions of 1 mm x 20 mm x 7.8 mm, 
where the ferroelectric c-axis is pointing along the 7.8 mm- 
long direction, and channels are aligned along the 20 mm- 
long propagation direction [Fig.[TJa)]. By choosing an ap- 
propriate thickness of the sputtered Ti stripes, and du- 
ration and temperature of in-diffusion, we form a sample 
possessing only a sing le allowed band [Fig. [jjb)-(d)]. This 
allows us to use standard end-facet coupling for study- 
ing linear and nonlinear light propagation in the lattice. 
Contrary to the single uniform lattice, the first band of 
our sample is split into two parts, separated by an addi- 
tional mini-gap that opens in the middle of the Brillouin 
zone. Finally, this sample possesses a defocusing saturable 
nonlinearity which is probably due to small impurity con- 
centrations incorporated into the substrate during high- 
temperature treatment, as well as intrinsic photorefractive 
defects like Nb on Li sites in congruently melting LiNbC>3. 
At slightly higher intensities or input powers [of the order 
of (10-100) /iW per channel] the observed nonlinear index 



changes of this sample are of approximately the same mag- 
nitude when compared to, for example, those fabricated 
on Fe-doped substrates. 
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Fig. 1. Geometry of the binary waveguide array (a), exper- 
imental output intensity distribution for homogeneous array 
excitation with Bloch momentum k = (b), output intensity 
distribution when a single element of the binary lattice is ex- 
cited (see arrows), either on the left (c) or the right hand side 

(d). 



To form dark solitons in our sample, we focus an ap- 
propriate light amplitude distribution onto the polished 
input facet, and record the temporal evolution of the out- 
coupled light by means of a 20 x microscope objective and 
a CCD camera connected to a computer. With the help of 
an additional plane reference wave we are able to monitor 
the phase distribution of the out-coupled light, too. For 
the case of a binary lattice different kinds of dark soliton 
solutions exist: a fundamental dark soliton in the semi- 
infinite gap, a bubble-like dark soliton that is located in 
the semi-infinite gap, too, and another dark soliton exist- 
ing in the extra mini-gap having an alternating phase dis- 
tribution. All solutions require excitation using a suitable 
amplitude profile and appropriate power to be formed. 
Because the nonlinearity of our sample is of saturable na- 
ture, we can, in a certain range, vary the input light power 
accordingly in order to generate the necessary nonlinear 
index changes for dark soliton formation. 

For the formation of a dark soliton in the semi-infinite 
gap above the center of the Brillouin zone at k = 0, we 
use a broad input beam with homogeneous intensity dis- 
tribution that is focused into the lattice with the help of 
a cylindrical lens. This beam passes a phase mask which 
imparts an additional phase of ir on half of the beam. The 
dark notch generated in this way is adjusted to coincide 
with the center of one element of the binary lattice. First, 
in Fig. [2|a) we show the output intensity distribution, 
when the lattice is excited on the input surface homoge- 
neously without using the phase mask. Adding the phase 
mask into the beam and using low input powers of the 
order of nW, we observe discrete diffraction of the gen- 
erated dark notch in Fig. |2|b). With increasing power, 
the sample starts to behave nonlinearly and the width of 
the dark notch decreases with increasing recording time. 
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After a few minutes it evolves finally into a dark soliton 
shown in Fig. [2|c), where the width of the dark notch 
matches the one on the input facet. This nonlinear state 
is robust and can be observed for more than one hour 
without noticeable changes. The corresponding phase re- 
lation is demonstrated in the interferogram presented in 
Fig. [2^d), which clearly shows the phase jump of it in the 
center of the structure. The generation of a dark soliton 
leads to a change of the corresponding refractive index 
profile and the formation of a positive defect located at 
the position of the dark notch. We test the guiding prop- 
erties of this defect by launching a broad probe beam of 
low intensity having a plane phase front, and observe trap- 
ping of light inside the recorded defect in Fig. [2^e) . The 
experimental setup used in this work is similar to the one 
in Ref. [12], but with the two-beam interference on the 
input side replaced by a phase mask. For the formation of 
a dark soliton in the semi-infinite gap above the center of 
the Brillouin zone at k = 0, we use a broad input beam 
(diameter 6 mm) of wavelength 532 nm with homogeneous 
intensity distribution that is focused into the lattice with 
the help of a cylindrical lens (/ = 100 mm). This beam 
passes a phase mask (half- wave MgF2 layer evaporated on 
glass substrate) which imparts an additional phase of tt 
on half of the beam. 

When the phase mask used to add the extra phase 
jump of 7r on the input light distribution is misadjusted 
(tilted), a dark notch (generated in part by a shadow re- 
gion induced by the mask) with both sides of the structure 
being in phase is formed. Using this input light distri- 
bution on the input facet, we are able to experimentally 
observe another nonlinear localized mode: The generated 
bubble solution has a similar intensity distribution as the 
dark soliton [Fig. |2jc)], but without a phase jump of tt in 
the center. The corresponding phase distribution on the 
output facet is given in Fig. [2](f ) . Note that in the ex- 
periment this nonlinear state is found to be unstable for 
long recording time, but may be observed during the ini- 
tial time of recording only. For longer recording or when 
the input power is further increased, a transformation to 
the dark soliton solution described previously occurs. This 
experimental observation suggests that the dark soliton 
should be structurally more robust than the bubble wave- 
form. 

To form a dark soliton in the mini-gap one needs to 
excite each dimer of the lattice in phase, with alternating 
phase jumps between neighboring units. Additionaly, a tt 
phase jump inside the central element is required to gen- 
erate a dark notch. By using an interference pattern with 
a modulation adjusted to fit to the grating period of our 
lattice, and by passing this beam through an additional 
phase mask for the generation of the dark notch, a suitable 
input light distribution has been generated. However, due 
to limited precision in the adjustment (this was mainly 
the limited symmetry in exciting each dimer with exactly 
equal powers in both channels), no clear discrete diffrac- 
tion pattern was observed for this situation on the output. 
Thus, we were unable to perform the corresponding non- 
linear experiments. A possible solution might be the use 



of a spatial light modulator to create a more precisely de- 
fined input profile, which was not available presently in 
the experimental setup used. It should be noted here that 
these "staggered" states have, in fact, been accessible to 
our numerical investigations. Nevertheless, the lack of ex- 
perimental control over their generation indicated above, 
as well as equally importantly, our generic observation 
of their large scale instability and the rapid evolution of 
their unstable background towards lattice dynamical tur- 
bulence have precluded us from considering these states 
further herein. 




Fig. 2. Light distribution on the output surface of the lat- 
tice by homogeneous excitation on input surface (a), discrete 
diffraction (b), dark soliton (c), interferogram of the dark soli- 
ton showing the tt phase jump (d), guiding properties of a dark 
soliton when the lattice is excited with a homogeneous probe 
beam (e), and interferogram of a bubble soliton when both 
sides of the dark notch are in phase (f). 



3 Theoretical Model and Analytical 
Considerations 

In order to theoretically model a setup corresponding to 
our experimental configuration, we consider a set of dis- 
crete nonlinear Schrodinger (DNLS) equations with a sat- 
urable onsite nonlinearity (see also [12 ) in the form: 

rlA I A -I 2 

i^ + ^+C^O + oJSj^B^O (2) 

The defocusing nature of the nonlinearity is encapsulated 
in the value of a = — 1; furthermore, we have used a value 
of n — 5 x 10 -4 which has been found to accurately ac- 
count for the saturation effect. The electric field of interest 
E n is defined from Aj and Bj as: 

E 2 j = Aj, E 2 j+i = Bj. (3) 
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Within this model, we seek stationary solutions in the 
form: 



A n = exp(if3z)a n , B n = exp(i/3z)6 n . 



(4) 



For the experimental waveguide array, by matching infor- 
mation about the linear diffraction properties of the lat- 
tice, we have inferred that C\ = 0.43 mm -1 and C 2 = 0.14 
mm -1 provide reasonable approximations to the experi- 
mental diffraction pattern observations [see Fig. l(c,d)]. In 
light of that, we have fixed the value of C\ to 0.43 mm -1 , 
but in order to appreciate the effect of variation of the 
spacing, we have looked at the relevant phenomenology 
as a function of the parameter C 2 - Furthermore, we have 
varied the propagation constant /?, in order to explore the 
family of nonlinear stationary, potentially observable so- 
lutions. 

We have focused, in particular, on two different types 
of solutions, namely dark [T51H4] and bubble [24] solitons, 
chiefly with a uniform background and a ni b n G M Vn. 

Upon identifying the relevant solutions through a fixed 
point computation, we have proceeded to examine their 
spectral stability properties by considering a linearization 
analysis. To that effect, small perturbations [of order 0(5), 
with < 5 <C 1] are introduced in the form 



A n (z : x) — e 



i(3z 



B n (z, x) = e^ z B n , + S(R n e^ z + S*e 



, (5) 



and the ensuing linearized equations are then solved to 
0(6), leading to the following eigenvalue problem: 



(6) 

for the eigenfrequency uo and the associated eigenvec- 
tor (P n , Q n , R n , S n ) T , where Li, L 2 , are the following 
operators: 



L 2 (x) v 



2|x| 



(1 + M 2 ) 2 ' 



(1 + M 2 ) 2 ' 



H n x n = dx n + C 2 x n±1 



to examine the linearization properties of such a back- 
ground. This can be done by using the following form: 



)exjp(ikn) 



b n = </>exp(z/cn) 



(8) 



with k being the wavenumber of the background. In our 
case, we have restricted our considerations to k = (un- 
staggered background) and k = tt (staggered background). 
Introducing the relation above into Eq. Q we get: 



P- 



1 



e(k) , 



(9) 



where e(k) = C\ + C 2 cos(fc). It is important to note here 
that in order that <\> G R, the condition that needs to be 
fulfilled reads: 



f3e[-e(k),-v-e(k)] , 



(10) 



which, when k = 0, will turn out (in our numerical results 
below) to coincide with the existence interval for dark soli- 
tons with unstaggered background. More generally, that 
condition constitutes a necessary one for the existence of 
dark / bubble solitons. In fact, as we will see below bub- 
ble solutions only exist in a sub-interval of the allowable 
parametric interval. 

The dispersion relation of the linear excitations corre- 
sponds to the continuous spectrum that will be identified 
in the linearization around our solitary wave solutions. 
This relation can be identified by decomposing the per- 
turbations as {P n , Q n , R ni S n } = {P,Q, R, S} exp(zgn) in 
Eq. Q and deriving the resulting condition: 

co 2 ± (q) = e(k)[e(k)-2h(k)]+g 2 (q)±2[h(k)-e(k)]g(q) (11) 
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g 2 (q) = C\ + C\ + 2CiC 2 cos( 9 ) . 



(12) 
(13) 



It is easy to check that Eq. ( 11 ) predicts the existence, 
when restricting to positive eigemrequencies, of two bands 
of real eigenvalues lj± when k = 0; besides, that equation 
also predicts for k = tt the existence of a band uo- with 
real eigenfrequencies and another band cj + with imaginary 
eigenfrequencies. This amounts to the modulational insta- 
bility of the background [25] which is detrimental for the 
survival of the staggered states mentioned in the previous 
section. 



and with 

v = - x 10" 3 = -2. (7) 

Hi 

As the relevant solutions (stationary dark solitons which 
possess a phase jump across the density dip, as well as bub- 
bles with no phase jump across the dip) live against the 
backdrop of a constant density background, it is relevant 



4 Numerical Observations 

We now turn to numerical observations, as a way of illus- 
trating the existence and potential stability of the local- 
ized states under consideration above. 

Figs. [3][4] show examples of the profile and the spectral 
stability planes of dark and bubble solitons with unstag- 
gered backgrounds. It can be seen that as the propagation 
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parameter is varied, the discrete dark soliton solution may 
go from a regime of spectral stability to one of weak os- 
cillatory instability and finally to one of stronger expo- 
nential instability. The former (oscillatory) instability is 
associated with complex eigenfrequencies, while the lat- 
ter (exponential) instability with purely imaginary ones. 
Similar features are observed for the bubbles presented in 
Fig. [4] A fundamental difference appears to be that for 
dark solitons an eigenfrequency pair associated with the 
wave first moves along the real frequency axis. This leads 
to resonances with the modes of the continuous spectrum, 
yielding oscillatory instabilities until a point of maximal 
excursion i.e., maximal Re(cj). Thereafter, the real part 
of this eigenfrequency starts decreasing and eventually 
crosses zero. Thus, it becomes unstable as an imaginary 
eigenfrequency causing an exponential instability. On the 
contrary, the relevant eigenfrequency pair first becomes 
imaginary for bubbles, leading to an immediate exponen- 
tial instability. Then, upon a maximal excursion, it returns 
to the origin, and subsequently moves along the real fre- 
quency axis. This gives rise to an oscillatory instability (al- 
though additional oscillatory instabilities may also exist). 
These features are presented for different values of C2 in 
Fig. [5] It is relevant to mention that in that figure in addi- 
tion to the experimentally relevant value of C2 = 0.14, pre- 
sented both for dark solitons (top panel) and for bubbles 
(bottom panel), cases of considerably smaller and larger 
C2 are examined in the middle panels. For smaller C2, the 
range of oscillatory instabilities shrinks (chiefly because 
of the narrower interval of continuous spectrum and the 
wider spectral gap in such a case). On the other hand, 
for larger C2, the corresponding instability is exponential 
and is found to arise throughout the examined interval of 
propagation constants (see third panel of Fig. [5|. More 
specifically, for C2 = 0.05, the exponential instability in 
the right panel of Fig. [5| only exists for the dark solitons 
with f3 > 0.81; for C 2 = 0.14, it exists for f3 > 0.48 and for 
C2 = 0.5, it exists for all j3 > —0.92. On the other hand, 
for the bubbles of the bottom right panel of the figure, the 
exponential instability is only present for (3 < 0.54 (notice 
the opposite sign of the inequality as discussed above). 

The detailed existence range of dark and bubble soli- 
tons is shown in Fig. [6j It is worth remarking that our 
numerics show that dark solitons always exist in the range 
f3 G [— (Ci + C2), — (Ci + C2 — v)\\ however, the existence 
range of bubbles is limited (to a narrower sub- interval), 
and these solutions cease to exist for C2 > 0.34. The 
stability properties together with the relevant instabil- 
ity growth rates are summarized in Fig. (T) It should be 
noted that there exist windows of stability inside the re- 
gion where the solitons are oscillatorily unstable; this fact 
is caused by the finite lattice size (see e.g., [9J for a rel- 
evant discussion) and is ubiquitous for dark solitons in 
lattices. Furthermore, this phenomenon disappears in the 
limit where the length of the lattice tends to infinity. 

In Figs. |8][9j we show typical examples of the dynami- 
cal evolution of the relevant waveforms, which have been 
obtained by a 4th order fixed-step Runge-Kutta method. 
In these examples, we typically append to the unstable so- 




Fig. 3. (Left) Profile of the normalized electric field and the 
corresponding intensity pattern for a dark soliton with C2 = 
0.14 and f3 = -0.54 (top), f3 = -0.27 (middle) and f3 = 1.13 
(bottom). The corresponding right panels show the spectral 
stability plane. In the latter, the existence of eigenfrequencies 
with non- vanishing imaginary part is tantamount to instability 
with a growth rate equal to the absolute value of this imaginary 
part. 



lutions a random perturbation of amplitude ~ 10 -3 , so as 
to seed the instability (through the projection of this ran- 
dom perturbation to the most unstable eigenmode of the 
solution). The right panel of those figures is restricted to 
twice the sample length, i.e., z = 40 mm (in order to gauge 
the relevance of the reported phenomenology to the exper- 
imental observations). Our general conclusions based on 
these dynamical runs (as well as others not shown herein) 
is that in most cases, the exponential instabilities may 
lead to a breathing of the light intensity (i.e., |£? n | 2 ) or to a 
potential destruction of the pertinent dynamical state. An 
example of the former type is shown in the bottom panel of 
Fig.[8]and is justified by the fact that for different values of 
the propagation constant f3 (and all other system param- 
eters being identical), there exist other members of this 
family of solutions which are, in fact, dynamically stable. 
An example of the latter type is shown in the top panel of 
Fig. [9] for the pertinent bubble state. On the other hand, 
instabilities associated with complex quartets of eigenfre- 
quencies are found to give rise to an oscillatory growth 
that is eventually seen to typically lead to propagation of 
excitations through the lattice with a non-zero angle (i.e., 
moving solitary waves). An example of this type can be 
seen e.g. for a dark soliton in the top panel of Fig.[5J while 
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Fig. 4. (Left) Profile of the normalized electric field and the 
corresponding intensity pattern for a bubble soliton with C2 = 
0.14 and ft = 0.10 (top panel), ft = 0.63 (middle panel) as well 
as /3 = 1.2 (bottom panel). The corresponding right panels 
show the spectral stability plane. 



for a bubble a similar observation but leading to a more 
complex dynamical state can be seen in the bottom panel 
of Fig. [9j Similar oscillatory instability evolutionary out- 
comes for a homogeneous lattice can be found e.g. in [TQl 

eng. 

As concerns the bubbles, we should note that their dy- 
namical instability in the cubic-quintic problem was sys- 
tematically studied in the works of [27,28 through a se- 
ries of analytical and numerical tools. There, it was found 
that stationary bubbles were unstable through a splitting 
towards stable traveling ones. The fundamental difference 
encountered here e.g. in Fig. [9] with respect to that setting 
is that in the present work, the state of vanishing inten- 
sity is dynamically unstable and hence the instability is 
more likely to break up the bubble into shallower propa- 
gating excitations, rather than to an expanding nucleus of 
vanishing intensity. This is what is observed through the 
asymmetric process (seeded by random noise) in Fig. [9J 
although it should also be added that in the presence of 
discreteness, the existence of the so-called Peierls-Nabarro 
barrier renders it unlikely that these excitations will main- 
tain their speed due to the emission of small amplitude 
radiative wavepackets, as discussed e.g. in [29] . 
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Fig. 5. Dependence with respect to /3 of the real (left) and 
imaginary (right) part of the linearization eigenvalues for some 
selected values of Ci- 



5 Conclusions 

In the present work, motivated by the experimental ex- 
amination of binary lattices with alternating couplings be- 
tween two distinct values C\ and C2, we have proceeded to 
theoretically model and subsequently analyze such states 
in the context of optical waveguide arrays. The two pro- 
totypical states under consideration have been discrete 
dark solitons and bubble solutions. The distinctive feature 
between these two waveforms is the existence (for dark- 
solitons) or non-existence (for bubbles) of a phase jump 
across the density dip associated with the solution. Fixing 
one of the two couplings, we have chosen to vary the other 
coupling, as well as the propagation constant of the solu- 
tion (as a measure of the nonlinearity) and to examine the 
existence and stability properties of such states. We have 
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Fig. 6. Existence range of dark solitons and bubble solitons. 
Dark (bubble) solitons exist in the region between the full red 
(dashed black) lines. 
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Fig. 7. (Top panels) Instability growth rates for dark (left) 
and bubble (right) solitons. The bottom panels show the type 
of instability experienced by the solitons, with the color cod- 
ing as follows: black — > stability; green [medium gray] — > oscil- 
latory instability; pink [light gray] — >> exponential instability; 
blue [dark gray] — >• oscillatory + exponential instabilities. 



found that dark solitons exist throughout the interval (of 
couplings and/or of propagation constant) in which the 
background state is found to exist. On the other hand, in 
the case of bubbles, the range of existence is found to be 
somewhat more limited, a feature that may be responsible 
for the experimental observation of their eventual "conver- 
sion" to dark soliton states. The existence of two tunable 
coupling parameters (or of the tunability of those in con- 
junction with the propagation constant) is found to offer 
some novel properties to the present model in connection 
to its standard homogeneous dynamical lattice counter- 
part. Firstly, the linear spectrum consists of two bands 
with a "mini-gap" between them, which is controlled by 
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Fig. 8. Propagation of the unstable waves of the middle and 
bottom panels of Fig. [3] respectively, when a random pertur- 
bation of amplitude ~ 10 -3 is introduced. The right panel is 
a zoom of the left panel. In the left panel, the resolution of 
the propagation is z = 0.1 mm, whereas in the right panel, the 
resolution is z = 0.01 mm. 
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Fig. 9. Propagation of the unstable soliton of the top and mid- 
dle panels of Fig. [4] when a random perturbation of amplitude 
~ 10 -3 is introduced. The right panel is a zoom of the left 
panel, as before. In the left panel, the resolution of the propa- 
gation is z = 0.1 mm, whereas in the right panel, the resolution 
is z = 0.01 mm. 



these parameters. This enables the existence of additional 
stationary (staggered) states whose dynamical instability 
has been briefly touched upon herein. Furthermore, the 
parameters enable transitions between oscillatory insta- 
bilities, dynamical stability and exponential instabilities 
for the same type of state, a wealth of possibilities that is 
not accessible in the standard DNLS model case. 

It would be particularly interesting to extend the present 
considerations to higher dimensional settings and espe- 
cially to two-dimensional arrays where different tunabili- 
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ties of the spacings may be accessible. On the one hand, 
there exists the tunability between different lattice direc- 
tions, essentially amounting to anisotropic nonlinear dy- 
namical lattices (see e.g. the early examinations thereof 
in [26], where interesting phenomena emerged from the 
breaking of the perfect symmetry of the square lattice). 
Yet on the other hand, there is the possibility of combin- 
ing this anisotropy with binary couplings, creating possi- 
bilities for inter- and intra-directional heterogeneity and 
their interplay. In that regard, the effect of such heteroge- 
neous nonlinear dynamical lattices on fundamental higher 
dimensional excitations such as discrete vortices would be 
particularly interesting to explore in future work. 
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